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ABSTRACT 

When an oscillating piston forces the gas in a closed 
tube to vibrate at a finite amplitude near an acoustic reso- 
nance frequency, shock waves are generated which travel back 
and forth in the tube. The behavior of these waves in tran- 
sition from the acoustic into the shock region for air at stan- 
dard temperature and pressure was studied at fundamental reso- 
nance frequencies up to 1000 cycles. An electrodynamic vibra- 
tor, capable of variable acceleration up to 140 meters per 
second per second, was used as a piston driver; it produced 
acoustic pressures in the tube up to 160 db sound pressure 
level. 

The dependence of pressure Wave shape on piston motion 
was investigated; it was discovered that over a large range of 
piston accelerations, nearly identical wave shapes were pro- 
duced at all frequencies provided that the same piston accele- 
ration was maintained. 

Although the largest shocks occur when the frequency is 
adjusted so that the shocks arrive at the piston just as the 
piston reaches its maximum withdrawn position, shocks can be 
observed at appreciable fractions of a period earlier or later. 
The effect of small changes in frequency away from resonance 
on this phase difference was studied, and a similarity para- 
meter was determined for this variation. 

Shock strength was measured as a function of both the 
phase of shock occurrence at the piston and piston accelera- 
tion. It was also found that acoustic resonance and shock 


resonance for a given tube occur at exactly the same frequency. 
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1. Introduction 
When an Oscillating piston forces the gas in a closed 
tube into vibration, a number of interesting effects occur, 
If one undertakes to study the motion of the gas in sucha 
tube according to acoustic theory, results are obtained which, 
although reasonably correct in most instances, predict particle 
velocities and pressures Which tend to infinity at resonance, 
Consider a cylindrical tube, closed at one end by a rigid 
wall and at the other end by a vibrating piston. Assume that 
the diameter of the tube is sufficiently small so that the 
waves travel down the tube with plane wave fronts. With particle 


velocity, u, as the dependent variable, the Wave equation is 





au _ 1 dtu 
Ces ef aes 
Whose solution may be written 
Ks t+kex) 
u(%t) = u,e!!% Rt) oy ae’ (1) 


Where ® = ke, and k = 2Tr/r » With © the angular frequency 


and A the wavelength, The boundary conditions may be approxi- 


mated by 
At x = 0, u(O,t) = wee, so that 
ee a es = Uo ) (2) 
ate. = | wUlbet)e=.0, so that 
—~ikL ikL 
Use ee ae =— © : (3) 
solving : 
a Sin kl =3 


If position is fixed at some particular value of x and if L is 
1 





held constant, then, as frequency is varied, both the numerator 
and denominator of (4) will vary. When kL is some multiple 

of TT, the particle velocity becomes infinite, except at points 
where k(L = x) is also a multiple of Tr, This leads to an ex- 


pression for the length of a closed tube at resonance: 
= a 4 ieee Se ee 


Ll e : 
It may be verified by experiment that even though the ampli- 





tude of oscillation of the driving piston is infinitesimal 
compared to the length of the tube, large-amplitude velocity 
and pressure vibrations are obtained when the length of the 
tube is as given above, It is to be noted, however, that the 
approximate expression given by (4) becomes indeterminate at 
the ends of a resonating tube. B. Miller and L. Olsen [12] 
make a more nearly exact analysis based on the wave equation 
written with a velocity-dependent damping term. They develop 


an expression for the length of a resonant tube given by 


U= so Be 52 7 = 


where & is the damping constant. In the usual experimental 
situation encountered in acoustics, the second term is negli- 
gible in comparison with the first, and the expression reduces 
to that given in the approximate analysis. 

In practice, of course, although particle velocities and 
pressures become very large at resonance, the effects of vis- 
cosity and heat dissipation are such as to limit the amplitudes 
to flnite values. If the piston amplitude is truly infinitesi- 
mal, the wave shape at resonance retains its sinusoidal forn, 


and the wave may be referred to as an acoustic wave, 





As the amplitude of piston oscillation is increased, 
the behavior of the gas changes. The most striking pheno- 
menon occurs when the frequency of vibration is near one of 
the acoustic resonance frequencies. This is the formation of 
shock waves which travel back and forth in the tube. With 
these finite amplitude vibrations, the pressure and velocity 
waves make a radical departure from the sinusoidal character 
observed at smaller amplitudes, and it becomes obvious that 
the acoustic theory represented by (4) no longer applies. 

Up to the present time, a complete theory describing 
Such finite amplitude vibrations has not been formulated. 
Several investigations of varying thoroughness have been con- 
ducted, which have led to somewhat inconclusive results, and 
it appears that more experimental data are needed. Frevious 
experiments were conducted using mechanically driven pistons 
which of necessity moved with constant amplitude, since with- 
out special linkages it is somewhat awkward to change the ampli- 
tude of such a driver. The experiments which are the subject 
of this thesis were conducted using an electrodynamic vibrator 
Which was easily variable over a large range of amplitude and 
frequency, so that measurements could be made under many dif- 
ferent conditions. This study was conducted toward a threefold 
purpose; 

a. To show the effects of variation of amnlitude on the 
ocCurrence and character of periodic shock waves. 

b. To increase our knowledge of the effects of other 


parameters, 








c, To demonstrate the feasibility of the use of an 
electrodynamic shaker for generating Shocks for further ex- 


periment. 


2. Historical Background 

Investigations into the phenomena of nonlinear gaseous 
vibrations in a closed tube were begun in 1935 by E. Schmidt 
it | and ©. Mayer-Schuchard iar who were interested in deter- 
mining the nature of this type of vibration as experienced in 
internal combustion engines. They used a tube (length 1200 cn, 
radius 3.4 cm, piston stroke 3.4 cm) filled with air at nor- 
mal density. A motorcycle engine was driven by an electric 
motor and its piston furnished excitation to the tube. Near 
the acoustic resonance frequencies, pressure-time histories 
measured at the closed end of the tube showed marked discon- 
tinuities. Between these pressure jumps the pressure records 
followed a curve resembling a sinusoid. Schmidt recognized the 
Jumps as shock waves. 

Using the same apparatus, E. Lettau |3| in 1939 made 
numerous observations of both pressure and velocity at several 
stations in the tube. His measurements were made by electro- 
mechanical means, and his records were obtained by photograph- 
ing the trace made by a moving light beam. He found that the 
Shape of the velocity and pressure profiles was strongly de- 
pendent on the location in the tube at which the measurement 
was made (see Fig. 1), and postulated an explanation involv- 
ing the summation of two partial waves traveling in opposite 
dlrections. 

Attempts to solve for the motion of the gas in a closed 
tube by mathematical means were begun somewhat later. J. Kel- 
ler lal, in 1952, obtained exact solutions for the one-dimen- 


éional gas dynamic equations, assuming an inviscid, non-heat- 








PRESSURE PARTIOLE PARTICLE 
VELOCITY DISPLACEMENT 


Figure 1. 


Variation over one cycle of pressure, 
particle velocity and particle displace- 
ment at different stations in a closed 
tube resonating at its fundamental fre- 
quency. (After Lettau [3] ) 





conducting medium with a particular equation of state, often 
referred to as a Chaplygin gas. One of the results obtained by 
Keller is that shocks should not occur in the tube unless the 
piston amplitude is in excess of A/2 , where A is the wave- 
length. That experiment produces shocks at piston displace- 
ment amplitudes of considerably less than this magnitude shows 
that in this respect, at least, Keller's solution does not 
accurately describe the phenomenon, 

R,. Betchov [6] made a theoretical analysis of the same 
problem, this time assuming a tube filled with a perfect gas, 
and taking into account the effects of viscosity and heat con- 
duction. He chose to make linear approximations of the motion 
of the gas, which he stated are valid in every portion of the 
tube except in the vicinity of the shocks themselves. A dis- 
continuity between the linear approximations represents the 
existence of shocks. From his solution he prepared charts of 
the variation of velocity and density which show qualitative 
agreement with the results of Lettau. 

R. Saenger and G. Hudson ale [8 | made a complete theore- 
tical analysis based on the idealized physical model of a gas 
element shown in Fig. 2. This shows the element together with 
the boundary layer of gas between it and the tube wall. The 
element moves as a unit under the action of pressure forces on 
its opposite faces and frictional retarding forces at the rim. 
The wall acts as a heat reservoir and the use of Newton's law 
of cooling permits the development of one-dimensional energy 
balance relationships, Saenger and Hudson assumed that gas 


motion at resonance may be represented by the sum of a finite 
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Figure 2. Idealized physical model of gas element 
used by Saenger and Hudson for their theo- 
retical development. 


continuous component and a finite discontinuous one, both of 
Which are periodic in time. These were obtained using the non- 
linear expressions derived from the basic flow equations, the 
Rankine=Hugoniot shock conditions, and the boundary conditions 
on the medium. They were then used to make a theoretical deter- 
mination of shock strength and mean temperature as functions 
of piston amplitude, tube length and dissipation parameters. 

Comparison between theory and the meager amount of experi- 
mental data available at the time showed qualitative agreement. 
One difficulty which faced Saenger and Hudson was the determina- 
tion of the proper dissipation parameters for heat conduction 
and frictional effects. They were unable to do this theoreti- 
cally, and were forced to use empirical values to obtain agree- 
ment between theory and experiment. 

Among the conclusions obtained in their treatment is the 
following: for piston amplitudes less than a critical value, 
periodic shocks cannot exist. Thus for a tube of fixed length 


driven at piston amplitudes below this critical value, it may 





be expected that a continuous, though not necessarily acoustic, 
treatment is appropriate. This implies that the transition from 
the acoustic region to the shock region should be a sudden pro- 
cess, with the shock suddenly appearing fully developed in place 


of a wave of acoustic proportions. 





3. Description of Periodic Shock Wave Phenomena 

In order to observe shock wave phenomena in a closed tube, 
it is necessary first to match the driver frequency with the 
length of the tube for acoustic resonance, This length is some 
integral multiple of A/2. For small values of piston dis- 
placement only the pure tone characteristic of acoustic vibra- 
tions is audible. But as the amplitude is increased, a buzzing 
sound is superimposed which indicates the presence of periodic 
shock waves. If the pressure-time history is observed on an 
oscilloscope, one notices the steep discontinuities occurring 
at the same time in each cycle. If the length of the tube is 
varied by a few millimeters, or the frequency by a few tenths 
of a cycle, the shape of the pressure curve changes drasti- 
cally, and the shock wave changes in magnitude. If either fre- 
quency or length are changed by more than minimum values, the 
shock condition disappears and the tube returns to an acoustic 
mode of vibration, which gradually diminishes in strength as 
the tube departs further from the resonant condition, 

Fig. 3 is an x = t diagram showing the motion of the shock 
wave as it travels in the tube, with respect to the motion of 
the piston, From this it can be seen that the shock wave reaches 
the piston at the instant that it begins the instroke portion 
of its cycle. In actuality, the frequency-tubelength combina- 
tion may be adjusted so that the arrival of the shock at the 
piston is considerably out of phase with respect to this maxi- 
mum outstroke position. It must be mentioned that the use of 
straight lines on the diagram of Fig. 3 is probably an ideali- 


zation, since it is likely that they should be curved, repre- 
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Figure 3. Space-time diagram for first resonance (f, = 
c/2L) showing the characteristic lines of slope 
= 1/c along which periodic shocks travel. 

senting changes in shock velocity with local conditions [6] : 

Particle velocity measurements can also be used to show 
the existence of shock waves, as can be seen from Fig. 1. 

The apparatus used in the present study proved to be very 
well suited to the study of periodic shock waves, primarily 
because it allows the precise setting of parameters and the 
quick determination of their effects on the shock waves. The 
results obtained can be more easily understood in terms of the 
following parameters and concepts, 

a. Resonance 

The concept of resonance is basically an acoustical one, 
indicating the condition at which the tube length is precisely 
adjusted to be nA/2, causing the sound waves in the tube to be 
of maximum strength. How the acoustic resonant condition is 
related to the production of maximum shock strengths had not 


been determined prior to the present study. 


11 


b. Piston Motion 

In the process of making theoretical developments concern- 
ing periodic shock waves, previous writers have expressed re- 
sults in terms of piston amplitude. Keller and Betchov use 
piston displacement amplitude in its dimensional form, while 
Saenger prefers to use a non-dimensionalized form which is essen- 
tially proportional to piston velocity amplitude. The results 
of Saenger are derived in terms of the quantity X/L, where X 
is piston displacement amplitude and L is the length of the 
tube, seeking to make theoretical predictions which are appli- 
cable to any length of tube, so long as it is resonant to the 


driving frequency. But 


& 
and for a resonant tube, 
pe ee ae 
ee 
Hence, 
ee ee 
ce TC Te (5) 


showing that the parameter X/L is actually proportional to pis- 
ton velocity amplitude. Experimentally, the measurement of dis- 
placement amplitude presents some difficulty, and it is often 
easier to measure piston velocity and integrate to obtain dis- 
placement amplitude, which is a simple matter if the piston mo- 
tion is sinusoidal. 

Considerable evidence was obtained from the present study 
that piston acceleration is a more appropriate parameter. Ac- 
celeration can also be calculated simply from piston velocity 


Measurements. 
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c. The parameter yD 

One of the most critical variables affecting the strength 
of the shocks and the shape of the waves is the phase difference 
between the time of the arrival of the shock at the piston and 
the time that the piston reaches its maximum outstroke position. 
This parameter, A: is illustrated in Fig. 4. In this illustra- 
tion the smooth curve is a velocity output from an electromag- 
netic pickup attached to the piston armature. When this curve 
crosses the horizontal axis, the piston is at one of its ex- 
treme positions. The phase of the extreme position can be mea- 
sured with respect to the occurrence of the shock. In the pre~ 
sent experiment, the measurement of velocity was made at the 
piston, and the measurement of pressure at the opposite end of 
the tube. One must assume that the shock takes half a period 
to travel the length of the tube in each direction. The para- 
meter A has been assigned a negative sign if it is such that 
the piston reaches its maximum outstroke position before the 
shock appears at the piston, and a positive sign if the piston 
extremum occurs after the shock. 

d. Tube shape and dimension. 

It seems likely that tube shape and dimension would have 
considerable effect on the behavior of the resonating’ air col- 
umn. The theory of Saenger and Hudson involves the effects of 
viscosity and heat conduction, both of which would be sensitive 
to changes in the tube dimension. If schlieren-optical tech- 
niques are to be used, it would be necessary to change the 
Shape of the tube to a rectangular cross-section. In ‘the pre- 


sent study, only one tube was used, and that of cireular shape. 


13 





"uoTtyTsod Yoous ayy 04 yZoedser urTM padns 

-eoau ‘Jutod (qyuemaoetdstp uUnuUtXxeu) AYFOOTSA OL8azZ ay4 Jo uot TSod ayy 

ST V dajeuered oyg “L[ytoo dnyotd AYyTooTea aya fo qndqyno 3u4 ST sarTnO 
yyoous suL °*sdo oog = IF 4e Usye. Bam suerZoTTToso asaul 


“OnTeA UF peysTys 
ST Y ce youeted out Se adeys aAem aanssedd fo uoty ete, ‘bh aIna Ty 


osve = V7 “0 po ay a 0(0S- = VY °e 








4, Description of Equipment 

ag Piston driver. 

For the study of threshold behavior of periodic shock waves, 
the desirable characteristics of a driver include (1) suffi- 
ciently high amplitude for the production of shocks, (2) con- 
tinuously variable frequency over a wide range, (3) easily ad- 
4usted and continuously variable piston amplitude, and (4) pure 
Sinusoidal motion. Although it is theoretically possible to 
obtain all of these qualities in a mechanical driver, design 
difficulties arise which would ultimately result in an expen- 
sive and cumbersome apparatus which would be difficult to use. 

The electrodynamic vibration shake table was decided upon 
as a driver closely approaching the four desirable states men- 
tioned above. For this study the equipment used was a Calidyne 
Model A88 Shaker, manufactured by the Calidyne Company, Winches- 
ter, Mass, This shaker is designed with the objective of pro- 
viding a light weight, yet rugged armature frame for transnit- 
ting vibratory forces to an attached load. A set of trunnions 
allows it to be tilted for horizontal operation. Tapped holes 
are provided for connecting the load to the armature of the 
Shaker. Into one of these holes a piston was screwed which was 
used to excite the air column in the tube. 

Included with the shake table were a Power Supply Model 68, 
an audio oscillator, a Power Amplifier Model 68, and a D. OC. 
Field Supply Model 102C, all Calidyne equipment. FREQUENCY 
and AMPLITUDE controls on the audio oscillator, as well as a 
GAIN control on the power amplifier, allowed variation of fre- 


quency from 5 cps to 10 kc, and of rms acceleration from zero 
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to about 140 meters per second per Beene: In addition, choice 
of proper impedance match (using IMPEDANCE MATCH selector on 
Amplifier Model 68) produced sinusoidal piston motion with less 
than one percent harmonic distortion in the velocity signal at 
frequencies up to 2000 cycles. Specifications of the Model A88 
shaker are included in Appendix I. 

The Model A88 Shaker is equipped with a velocity pickup 
which consists of a coll attached to the table of the shaker 
so as to move in the field of a permanent magnet attached to 
the shaker body. The voltage produced is proportional to the 
relative velocity between the shaker table and body. This vol- 
tage can be continuously monitored as an indication of the ab- 
solute table velocity. The use of ana. c. voltmeter allows 
the shaker to be adjusted for operation at desired velocity am- 
plitudes. 

b. Microphone 

An Altec 21-BR-200 microphone system with Amplifier Type 
526B was used as a pressure pickup. This microphone was used 
because of its rugged construction and essentially flat frequen- 
cy response curve. Its relatively aioe rise time of about 60 
microseconds placed a limitation upon the accuracy of observa- 
tions made at higher frequencies. A comparison calibration of 
the microphone pressure sensitivity was performed using a Wes- 
tern Electric microphone Model WE-640-AA as a standard. The re- 
sults of this calibration are contained in Appendix II. 

ec. Tube and Piston. 

The study was performed using a 1-3/4 inch diameter alumi- 


num tube approximately seven feet long. A simple aluminum pis- 
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ton of approximately the same diameter was machined, and was 
connected by a short threaded rod to the shaker. Initially, a 
metal bracket was constructed into which the tube was to be 
threaded; the bracket was to be bolted to the face of the shaker 
with the intention of holding the end of the tube rigidly in 
place. Because of alignment difficwlties, trial runs were made 
with the tube loosely supported on a wooden frame, with the end 
floating on the vibrating piston. Using this method of support, 
it was found that the tube became almost self-aligning, and no 
other difficulties were immediately presented. Therefore the 
rigid support was discarded, and the method of non-rigid sup- 
port was used during the remainder of the investigation. 

It was also found destrable to make some changes in the 
piston as the experiment progressed. As originally constructed, 
the piston was smooth-walled and made metal-to-metal contact 
with the walls of the tube. Provided that the piston was kept 
well-lubricated (a mixture of Molykote lubricant and instrument 
o11 was used), this construction seemed to be fairly satisfac- 
tory, and was employed for several weeks. Using this piston, 
it was noted that changes in alignment occasionally caused un- 
due vibration and had some effect on the experimental results. 
On one occasion a combination of misalignment and insufficient 
lubrication led to extensive galling of the piston and tube, 
and it was realized that a redesigned piston was necessary. 

The damaged piston was machined down a few thousandths of 
an inch to remove the galling, and incidently reducing its dia- 
meter so that metal-to-metal contact would be avoided. A groove 


was cut into which an O-ring could be fitted, and the piston was 


17 





reinstalled on the shaker. The reconstructed piston proved to 
be more satisfactory than the old one, since it was even easier 
to align, and was more nearly gas-tight than before. The O- 
ring was Well lubricated with vacuum grease to keep wear to a 
minimum. 

The end of the tube opposite from the piston was termina- 
ted with a movable aluminum plug which was bored so that the 
microphone could be mounted flush with its face. Two O-rings 
provided an effective seal, and a long rod threaded into the 
rear of the plug allowed the length of the tube to be varied. 

d. Oscilloscope 

A Tektronix Model 545 oscilloscope was used as a recor- 
ding device because it possesses certain highly desirable 
characteristics. Since this study involved the measurement of 
steeply-rising pressure gradients, a fast oscilloscope rise 
time was essentials; the Tektronix scope has a 0.01 microsecond 
rise time. The use of the plug-in preamplifier unit type QA, 
Which provides a dual trace, allowed simultaneous observation 
of the pressure signal from the microphone and the shaker ve- 
locity signal. The delay feature of the Tektronix together 
With its sweep magnifier allowed measurement of phase differ- 
ences accurate to about two degrees. 

e. Frequency counter 

Because of the extreme sensitivity of the experiment to 
frequency, it was found necessary to measure frequency to the 
nearest tenth of a cycle per second. For this purpose, a Hew- 
lett-Packard Frequency Counter Model 5210 was used. Using a ten- 


Second count, the required precision in frequency was obtained. 
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5. Experiment and Results 

a. Preliminary Observations 

Once the equipment was set up, a series of preliminary 
observations over a wide range of amplitude and frequency were 
made. The following are among the facts immediately noticed: 

The transition from the acoustic domain to the shock do- 
main as amplitude is increased is a gradual one. Once the fre- 
Quency-tubelength combination is properly adjusted for the oc- 
currence of shocks, increasing piston motion from zero ampli- 
tude produces a sinusoidal pressure output which gradually dis- 
torts. The first indication of an incipient shock wave is the 
occurrence of a discontinuity in the first derivative of the 
pressure-time curve. Further increase in amplitude increases 
the bend until it becomes vertical. Still further increase 
brings about a break in the curve which is the characteristic 
appearance of the wave containing the shock. 

Although the shape of the pressure-time wave changes dras- 
tically during the transition from the infinitesimal- to the 
finite-amplitude region, wave shapes of almost identical ap- 
pearance can be produced at each frequency for a given A\ -value 
by selecting the proper piston amplitude. In addition, when 
the shaker is oscillating at its maximum available amplitude, 
the pressure wave shapes obtained at different resonant frequen- 
cles are approximately the same. Using the Altec microphone, 
the above-mentioned effect is frequency-limited, because as 
frequency is increased, the rise time of the instrument assumes 
a more and more significant fraction of the period of vibration, 


masking whatever similarity may exist. 
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The value of Aisa very sensitive function of both f and 
L. For a system adjusted so that A= 0, at a frequency of, 
say, 200 cycles, a change of two-tenths of a cycle will cause 
QV to change by about ten degrees, and a change of three cycles 
per second will drive AX beyond 90°, whereupon the shock disap- 
pears. Small changes in the length of the tube have a corres- 
ponding effect on LS 

Shocks are found to exist in the tube under the most fa- 
vorable conditions for their production (1. e., large ampli- 
tude driver motion) through a Z\-range of from about -110° to 
+909. As A passes outside of this range, the shock diminishes 
in size and finally disappears, although the vibration does not 
immediately become sinusoidal. 9 

b. Analysis of Wave Shape. 

The fact that similar wave shapes can be produced at each 
frequency by selecting the proper amplitude led to the supposi- 
tion that there may exist a parameter, independent of frequen- 
cy, by which wave shape may be specified. In order to pursue 
this supposition, a harmonic analysis of the wave shape was 
conducted for several values of amplitude and frequency. This 
analysis was made using a Hewlett-Packard Wave Analyzer Model 
302A. Several sets of data were taken for frequencies of from 
100 to 500 cps and for the second through the twelfth harmonic, 
as functions of rms piston velocity which was read in milli- 
volts output of the velocity pickup. These data were then 
plotted against piston displacement, velocity and acceleration, 
It was soon noted that, although the curves of percent harmonic 


versus displacement and versus velocity bore little similarity 
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to each other, the curves of percent harmonic versus piston 
acceleration were strikingly similar, and in some cases iden- 
tical, for all of the frequencies tested. It could then be 
postulated that the similarity parameter connecting the shocks 
is piston acceleration. 

In an effort to substantiate this postulate, a series of 
twenty-seven oscilloscope photographs were taken of the pres- 
sure-time wave as seen at all combinations of three different 
frequencies, three different values of Le and three different 
acceleration amplitudes. In order to make a visual comparison, 
the oscilloscope vertical sensitivity and sweep time were ad- 
Justed so that the wave was exactly the same height and the 
same length in all cases. 

Extraction of the data from the oscilloscope photographs 
was accomplished using a projection technique. Each photo 
was placed in an opaque projector and the curve projected onto 
a large sheet of graph paper. Using arbitrary numerical scales 
along the vertical and horizontal axes, the ordinates of the 
pressure-time curves were read and recorded for a number of 
equispaced points along the abscissa. Initially, seventy-two 
points were read off each curve, but this was later reduced to 
foc y 

The harmonic analysis itself was done using the CDC Model 
1604 digital computer installed at the U. S. Naval Postgraduate 
School, A program was written which would accept input data 
in the form of the data points from the curves, perform the 
Fourier analysis, and then divide each amplitude by the ampli- 


tude of the fundamental. In this way, the results were in 
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effect normalized to a common reference, in the same way as 
for the work done with the electronic wave analyzer. | 

The data obtained in this way were then plotted, with per- 
cent of the fundamental amplitude as ordinate, and root mean 
square acceleration, A as abscissa, If acceleration is indeed, 
as postulated, a similarity parameter, then for a given value 
of Ne the points taken for common A Should be the same re- 
gardless of frequency. The results obtained in this portion 
of the experiment are plotted in Figs. (6), (7), and (8), and 
show that agreement was obtained to within at least two per- 
cent in every case, and thus tend to confirm the basic postu- 
late. It is believed that departures from the anticipated agree- 
ment are due in part to small fluctuations in the value of A 
for the points of different frequency and to small changes in 
microphone sensitivity to the different harmonic frequencies, 
Even though perfect correlation was not obtained, it must be 
realized that 1f the harmonic amplitudes were plotted with re- 
spect to a parameter other than acceleration (such as velocity 
amplitude or displacement amplitude), the curves would be so 
dissimilar that they could not be plotted on the same piece of 
paper. 

In order to demonstrate further the meaning and effect of 
acceleration as a similarity paraneter, two sets of photographs 
taken at different piston frequencies are included herein for 
comparison. The first (Fig. 9) shows a pair of oscillograms 
of the pressure-time wave taken at nee cade value of piston 
velocity amplitude, and the second (Fig. 10) shows a similar 


pair taken at the same value of piston acceleration amplitude, 
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Figure 9. Oscillograms of pressure-time Waves at 
different frequencies with the same piston 
velocity. 
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Figure 10. Oscillograms of pressure-time waves at 
different frequencies with the same piston 
acceleration. 





The reader will note the remarkable similarity in the shape 

in the second case. The slight difference in height of the 
shockfront in the latter set may be due to the fact that micro- 
phone rise time assumes a relatively larger portion of the 
cycle in the higher frequency case, The small fluctuation 
superimposed on the basic waveform is due to the ringing of the 
microphone at the frequency of its mechanical resonance, 

It must be pointed out that merely determining the har- 
monic amplitudes of the waves does not in itself give a conm- 
plete picture of their construction, since the complete Fourier 
analysis is given by 

P= A, + Aycos(wt + g,) + A,cos(2ut + Bo) ts « 

and no determination of the constants Bs was made in this ex- 
periment. Indeed, the fact that the visual appearance of the 
wave changes drastically with A while the curves of Megs. § 
through 8 change only slightly at larger amplitudes indicates 
that there must be a substantial change in the g, as Zita 
changed. The difficulty in determining the g, is due to the 
-asymmetry of the waves which prevents the selection of a satis- 
factory zero reference which will be common to all waves, re- 
gardless of A. This same difficulty made it impossible to de- 
termine reasonable results for Ao; the computer program as 
written solved for this quantity, but since no common vertical 
zero reference was avallable for all waves analyzed, the numbers 
Obtained were meaningless. 

c. Variation of AQ with acceleration. 

In the course of the investigation it was noticed that 


during a series of measurements which were intended to be at 
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constant ZA, a change in the value of Awas often accompanied 
by a change in the value of AA, even though the frequency was 
maintained constant. Therefore, it was decided to investigate 
this variation. A series of measurements Was made which began 
with an initial setting of 2 at the maximum value of acceler- 
ation available on the shaker; A was then recorded as the ace 
celeration was reduced. This procedure was carried out for 
seven different initial A -settings for each of three differ- 
ent frequencies. The value of Zs measured at the maximum avail- 
able value of A will be termed Yar For each frequency sa 
Similar set of curves were obtained, included herein as Figs. 


11 through 13. In general, if the aS was positive, the 


max 
value of Z\ became more vositive as “4 was reduced; sSlmilarly, 
a negative eee, became more negative. However, it was inter- 
esting to note that eis zero did not remain zero with 
decreasing A », but rather became positive. In addition, there 


was one negative value of A which was different for each 


max? 
frequency checked, for which A\ remained constant as Al was re- 
duced. The measured positions of these constant- ZA lines are 
shown as dashed horizontal lines in Figs. 11 through 13. 

The close similarity of these sets of curves for the dif- 
ferent frequencies again demonstrates the existence of A as a 
Slmilarity parameter. 

Ge sNOCk Strength 

The most striking feature of the nonlinear wave of finite 
amplitude is the appearance of a sharp pressure discontinuity 


as the shock passes the recording device. When the microphone 


is installed in the rigid tube-end as it was in this experiment, 
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a double-amplitude shock is observed, as the incident and re- 
flected shocks are superimposed. 

Shock strength variation was studied over nine different 
values of Z\ at each of six frequencies, while varying A over 
its available range. The strength of the snock was measured by 
recording the height of the shock as observed on the calibrated 
oscilloscope face, From the vressure calibration of the micro- 
phone it was possible to obtain shock strength as a function 
of f, AX, and A. A series of quasi-linear curves were ob- 
tained which are included herein as Figs. 14 through 17. 

Although the response of the microphone 1s not perfectly 
constant over all frequencies, it is sufficiently flat over many 
harmonics that the single figure of -85 db re 1 volt per dyne 
per square centimeter can be used to compute pressure from volte 
age output. It was also necessary, of course, to halve the 
pressure obtained because of the double-amplitude nature of 
the wave as measured, in order to compute incident shock over- 
pressure. 

The qualitative characteristics of the shock variation are 
Similar for all frequencies studied, At any A -setting, as 
the acceleration is reduced, the shock diminishes in strength 
and finally disappears. At the time that the shock disappears, 
the wave retains the form of a distorted sine wave, and it is 
moe until A is much further reduced that the pure sinusoidal 
character of the acoustic wave is assumed. The value of for 
Which the shock disappears is not the same for all frequencies, 
being about twice as great for 1000 cycles as for 100 cycles. 


In order to test the previous conclusion that the shape of 


30 





the wave can be predicted from piston acceleration, each value 
of shock strength data was divided by the amplitude of the fun- 
damental component of the wave measured under the same condi- 
tions. This was done to take into account the varying strengths 
encountered at the different frequencies. The points resulting 
from this calculation were then plotted against A, giving the 
set of curves shown in Figs. 18, 19 and 20. These curves essen- 
tially support the contention that Al is a Similarity parameter, 
since for every A, the points obtained by the above operation 
fall, within the limits of experimental precision, along a sin- 
gle line, regardless of frequency. 

By a careful study of the work of Saenger and Hudson is], 
one may conclude that the existence of the similarity parameter 
A for shock strength can be inferred at least qualitatively 
from their predictions. Their development results in the ex- 


pression 


fees OMe aE 
R m7 L ABE lA 6) 


where © is the dissipation parameter for viscous losses and 
Bis a parameter derived from the proportionality constant in 
Newton's law of cooling. According to Saenger and Hudson, 
agreement between experiment and theory could not be obtained 
unless empirical values for B and ® were used. These were 


such that ®)7B. Hence (6) becomes, approximately 


AP X GY | 
NS ce, AO ee Sin CA = 
es Tels £ in ( ) (7) 


In a theoretical investigation of the propagation of sound 
in a closed tube, Rayleigh [9] derived an expression for the 


tangential force per unit area on the wall by viscous effects 
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of air moving at unit velocity. Expressed per unit length of 
tube, this is 


To obtain the ¢ of Saenger and Hudson from this, F must be di- 
vided by the mass of the gas per unit length, pra” (because of 
their use of the Lagrangian notation), and by a (because of 
their use of nondimensionalized quantities), giving 
$= et 
- C4 

Numerically, this does not agree with the empirical results of 
Saenger and Hudson. Nevertheless, the important consideration 
here is the frequency dependence; it is noted that Pf ls pro- 


portional to f72, Hence, a further approximation of (7) is 
2 x . 
RQ ae. L 
Applying (5), equation (8) may be written 


ra = cal NG sin (A-) 








Which becomes 


= = 1a Al sin (A~ £) (9) 


The behavior of P, 1s not well-known, either as a func- 
tion of A or of frequency. More investigation must be done to 
reconcile if possible the facts that, although AP is near maxi- 
mum at A= 0 (see Fig. 21), the right member of (9) is near 
minimum at the same value of A. 

e. Phase Diagrams 

In taking the data for the shock strength variation, the 
minimum value of A for which a shock could be measured was re- 


corded at each of the resonant frequencies studied. This mini- 
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mum occurred at or near the A=0 frequency for the tube. 
Plotted against frequency in Fig. 24, these points furnish a 
phase diagram which separates the A-t plane into regions in 
which shocks occur and do not occur. These are labelled "Shock 
Region" and "Acoustic Region," respectively. A similar phase 
dlagram, prepared for displacement amplitudes instead of ac- 
celeration, is included as Fig. 25. 

ir A Were in fact a perfect similarity parameter, then 
the curve of Fig. 24 would be a horizontal straight line. 
since it is not, this shows that although the concept of UE, 
as a similarity parameter for shock waves is valid over a large 
range of piston amplitudes, it breaks down in the limit of very 
small amplitudes. 

The theory of Saenger and Hudson predicts the minimum 
piston amplitudes at which shocks may occur. A calculation 
based on this theory for f,; = 100 cycles per second predicts a 
minimum displacement amplitude for the shock region which is 
about thirty times that obtained for the 100 cycle point on 
mie. 25. 

f. Check of Correlation Between "Acoustic" and "Shock" 
Resonance Frequencies. 

Although it has long been known that the production of 
periodic shock waves in a closed tube occurs at frequencies 
which are very near the acoustic resonant frequencies of the 
tube, the writer knows of no instance where it has previously 
been verified that these frequencies are in fact the same. In- 
deed, the difference in magnitude of X for the two cases, as 


well as the more complicated theory applicable in the finite 
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amplitude case, leads one to expect that they might conceivably 
be slightly different. 

A check on these frequencies was a simple matter with the 
variable-amplitude driver used. First, the frequency-tubelength 
combination was adjusted so that Z\ was exactly zero with the 
shaker operating at maximum acceleration. Then A was reduced 
using the GAIN control on the shaxer until the output of the 
microphone was near its minimum detectable value. The sinus-~ 
Oidal appearance of the microphone output assured that this 
level was well within the acoustic domain. During the process 
of reducing the gain control, the frequency counter indicated 
that frequency was maintained constant to the nearest tenth of 
a cycle per #econd. 

Finally, the microphone output Was measured on ana. ¢c. 
voltmeter. Several measurements were made following this pro- 
cedure, at frequencies of 100, 300, 500 and 100 cycles. In 
every case, the voltmeter was found to be sharply peaked when 
the final measurement was made, and the slightest change in 
frequency in either direction caused an abrupt drop in the 
sound pressure aS measured by the microphone, 

From this it is concluded that the acoustic resonance fre- 


quency of a tube and the A = O condition for periodic shock 


ma xX 
waves occur at exactly the same frequency, Within the limits 
of exverimental precision. 

(In view of the shift in the Amax = O lines of Figs. 11, 
12 and 13 as Al is decreased, the results are especially in- 
teresting. Examining the curves, one might suspect that the 


values of Lee for which ae - O, represented by the dashed 


34 





lines, rather than the Ano, = 0 condition would be found to be 
the determining factor for the acoustic resonance frequencies, 
but this is definitely not the case.) 

g. Variation of Aj, with frequency. 

During several series of measurements, the variation of 
La With frequency was studied. As could be anticipated, it 
was found that the change in frequency necessary to produce a 
given change in A,., was greater as frequency was increased. 
The frequency at which the tube is resonant to the fundamental 
Will be denoted by f,; then the change necessary to produce a 
given Qiox 15 expressed by the quantity f - fj. When f - fj 
was plotted against O,,,, a family of monotonically decreasing 
curves was obtained (Fig. 22). Since all curves are of similar 
Shape, a parameter was sought by which the variation with re- 
spect to frequency could be quantitatively determined. After 
a short exploratory process, it was discovered that the quan- 
tity (f - £,)/f;2 is the same function of Anax, for all fj. 
This result is shown in Fig. 235. 

h,. Thickness of Shock Wave. 

All oscillograms taken show the shock Wave as an abrupt 
rise in pressure which takes place in a time so short that it 
could not be measured by the Altec microphone with its rise 
time of about sixty microseconds. In an effort to determine 
more about the structure of the shock, a new microphone was 
constructed. The pressure-sensitive element of the new micro- 
phone was a barium titanate wafer with a two-megacycle resonant 
frequency. The wafer was mounted in a brass plug fitted in the 


same manner as the plug which held the Altec microphone, 
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Using a high-gain preamplifier with the Tektronix oscil- 
loscope, the output of the barium titanate crystal in response 
to periodic shock waves could be observed. Although there was 
a considerable amount of ringing (due to the mechanical reso- 
nance of the crystal) superimposed on the familiar profile of 
the shock wave, the pressure jump of the shock wave was clear- 
ly visible, and its rise time could be measured. Careful mea- 
surements at several frequencies using the expanded scale of 
the oscilloscope gave rise times of the order of one-half mic- 
mesecoOnd, From this it 1s concluded that the times measured 
were those of the instrument and not of the shock wave, and 
that the shock pressure jump must be completed in less than 
one-half microsecond. The shock thickness, then, must be 
less than 0.2 millimeter, which is about the distance a sound 


wave travels in one-half microsecond, 
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Figure 11. Variation of A as A. is reduced. f, = 100 cps. 
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Figure 18. Normalized Shock Strength versus A. 


Vertical scale is AP in volts (as measured on 
oscilloscope) divided by fundamental of wave in rms volts. 
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Figure 19. Normalized Shock Strength versus A ; 


Vertical scale is AP in volts (as measured on 
oscilloscope) divided by fundamental of wave in rms volts. 
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Figure 24. Phase Diagram Showing Threshold Piston 
RMS Acceleration for Production of 
Periodic Shocks 
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6. Conclusions 

The existence of acceleration as a similarity parameter 
has been definitely demonstrated. This fact is important when 
considering the use of an electrodynamic shaker as a driver for 
generating shock waves. such shakers are rated at a certain 
maximum force output, which implies a maximum acceleration. 
Since acceleration is a similarity parameter, this implies that 
if a shaker is satisfactory for generating shock waves at one 
frequency, it can be used for this purpose at other frequencies 
Within its available range. Because of the frequency limita- 
tion of the Altec microphone, the upper limit at which shocks 
can be generated was not found. 

An examination of the plots of harmonic analysis, A, and 
normalized shock strength with respect to A shows a tendency 
for these curves to level off at the higher accelerations. 

This fact can be used to indicate the limits of the threshold 
region for periodic shock waves. Unless unforseen excursions 
of these curves take place at higher values of A, it can be 
concluded that at accelerations of approximately twice that 
available with the equipment used in this experiment, the 
CUsveswoOnerieas 6.9). 6, Ieee ios. 19 and 20 will have 
flattened out. Since the Calidyne Shaker Model A83 has avail- 
able accelerations of the order of 20g, this means that a ma- 
chine with rated accelerations of the order of 40g - 50g will 


be of sufficient capacity for the study of threshold phenomena. 
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7. Suggestions for Further Experiment 

saenger and Hudson [3] State that a small amount of ex- 
perimental data shows qualitative agreement with their theory. 
Because the behavior of P, was not measured in the present ex- 
periment, no additional evaluation of the correctness of the 
theory can be made here, A further step in the present process 
would be YO-obtain= the variation of P, with ae and rae along 
with shock strength information, and use this information to 
obtain a quantitative evaluation of their theory. 

In this process it would be necessary to determine ex- 
perimental values for the constants B and ®, It is known 
that the experimental and theoretical values for B are in fair- 
ly close agreement, and that this is not true for f. One 
possibility would be to use the digital computer to compare 
experimental data with that calculated using various values of 

2, untila Pof optimum size to sult experimental fact could 
be found. Once this optimum f is found, it could be deter- 
mined whether the theory adequately describes periodic shock 
waves. 

The relatively slow rise time of the microphone employed 
in the present experiment limited the extent to which frequency 
varitation could be investigated. The microphone rise time be- 
came a Significant fraction of the period when frequency was 
increased to 1000 cycles per second, and no investigations 
above that frequency were considered worthwhile. It 1s recom- 
mended that an investigation be made of the effects encountered 
at higher frequencies, using a more suitable means of pressure 


measurement. If a faster-acting microphone 1s used to measure 


Se 





pressure, it will be necessary to use a suitable amplifier, 
with a filter to remove the mechanical ringing inherent in such 
equipment. 

Schlieren-optical techniques offer another means of inves- 
tigating shock phenomena oa. eee Precise measurements of 
density variations in the air column can be made, and, pro- 
vided that one additional parameter, such as particle velocity, 
is measured, the corresponding oressures can be obtained using 


the equation of state. 
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APPENDIX I 


SPECIFICATIONS OF MODEL A88 SHAKER 


Characteristic 
Rated Force Output 
Field Current Recuired 
Continuous Armature Current 
Force/Current 
Longitudinal Resonance 
Armature Weight 
otroke 


Velocity Signal Calibration 
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100 lbs @ 5.26 amp rms 
5.0 amperes 

5.3 amperes 

19 lbs/ampere rms 

2560 cps 

2.68 lbs 

t inch peak-to-peak 


0.110 vyolt/inch per Sec. 
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Microphone Calibration Curve 
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APPENDIX III 
TABLE OF SYMBOLS 
Root mean square piston acceleration 
Radius of tube 


Dimensionless form of proportionality constant 
in Newton's law of cooling 


Velocity of sound 
Driving frequency of piston 
Fundamental resonance frequency of tube 


Tube length measured from mean position of 
piston to closed end of tube 


Shock strength -- the positive jump in pressure 
across the shock, 


Instantaneous pressure at a point in the tube 
Where gas is at rest 


Gas particle velocity 

Piston velocity amplitude 

Coordinate of position along length of tube 
Piston displacement amplitude 

Ratio of specific heats, ¢cy/cy 

Phase lag of the piston in its outstroke (maxi- 
mum withdrawn) position, behind the appearance 
or the shock at the piston. 

The value of Q measured with the shaker opera- 
ting at its maximum available A. 

Wavelength 

Viscosity of gas 

Kinematic viscosity 

Density of gas 


Dimensionless form of viscosity parameter 


— [lv-1)B*—2(7-1)Bs - 41/[v- 8 + I 
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